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ABSTRACT 

Magnetohydrodynamic (MHD) kink waves are ubiquitously observed in the solar atmosphere. The propa¬ 
gation and damping of these waves may play relevant roles for the transport and dissipation of energy in the 
solar atmospheric medium. However, in the atmospheric plasma dissipation of transverse MHD wave energy 
by viscosity or resistivity needs very small spatial scales to be efficient. Here, we theoretically investigate the 
generation of small scales in nonuniform solar magnetic flux tubes due to phase mixing of MHD kink waves. 
We go beyond the usual approach based on the existence of a global quasi-mode that is damped in time due 
to resonant absorption. Instead, we use a modal expansion to express the MHD kink wave as a superposition 
of Alfven continuum modes that are phase mixed as time evolves. The comparison of the two techniques ev¬ 
idences that the modal analysis is more physically transparent and describes both the damping of global kink 
motions and the building up of small scales due to phase mixing. In addition, we discuss that the processes 
of resonant absorption and phase mixing are intimately linked. They represent two aspects of the same un¬ 
derlying physical mechanism: the energy cascade from large scales to small scales due to naturally occurring 
plasma and/or magnetic field inhomogeneities. This process may provide the necessary scenario for efficient 
dissipation of transverse MHD wave energy in the solar atmospheric plasma. 

Subject headings: Sun: oscillations — Sun: atmosphere — Sun: magnetic fields — waves — Magnetohydro¬ 
dynamics (MHD) 


1. INTRODUCTION 


Recent high-resolution observations indicate that 
verse magnetohvdrodvnamic (MHD) waves are ubicj 

trans- 

uitous 

throughout the solar atmosphere (e.g. .IDe Pontieu et alJ 

200 ^ 

Tomczvk et al.ll2007l: iLin et alJl2009l:IOkamoto & De Pontieu 

201 ll McIntosh et alJ 12011: Kuridzeetal. 

2012 ; 

De Pontieu et al. 120121 iMorton & McLaughlinI 

2013; 


Morton et al.l l2014[) . From the theoretical point of view. 


the observed waves are interprete d as MHD kink waves 
on magnetic flux tub es (see, e.g., lEdwin & RobertsI 119831 : 
iGoossens et akl |2012|) . Both observations and theoretical 
aspects of these waves have be en re cently reviewed by 
iDe Moortel & NakariakovI (1201 2^ and iMathioudakis et al] 
(I2OI3I) . It is believed that the propagation of transverse 
MHD waves and their dissipation may play relevant roles 
for the heating a nd energy transport in the solar atmospheric 
plasrn a ('see, e.g. |Cargili & de Moortelll201 U iMcIntosh et alJ 
|2^ iParnell & De Moortell 120121: iHahn & SavinI |2014t) . 
Heating may be caused by wave energy dissipation due 
to Ohmic diffusion and/or shear viscosity. However, in 
the solar atmosphere Ohmic and viscous dissipation are 
inefficient unless very shor t spatial scales are involved (see 
iGoedbloed & Poedti |2004) . For instance, the magnetic 
Reynolds number, R,„, indicates the efficiency of Ohmic 
diffusion. It can be defined as 


R 


m 


Iqvq 

n 


( 1 ) 


where /q and vq are characteristic length and velocity scales of 
the plasma, respectively, and 77 is the Ohmic diffiisivity. We 
can relate vo with the plasma Alfven velocity (~ 1000 km s“* 
in the solar corona) and Iq with the length scale of the wave 
perturbations. The Ohmic diffiisivity depends on the plasma 
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properties and is extremely sm all for atmospheric plasma con¬ 
ditions (see, e.g.. lPriesil2014l) . This results in R^ ~ 10'^ for 
typical parameters in the solar corona. However, /?„ < 1 is 
needed for efficient dissipation. Such small values of the mag¬ 
netic Reynolds number can only be achieved in the solar at¬ 
mosphere if Iq is very small. Thus, some physical process able 
to transfer wave energy from large spatial scales to small spa¬ 
tial scales is needed for Ohmic dissipation of transverse MHD 
wave energy to be efficient in the solar atmospheric plasma. 
The same result holds in the case of dissipation due to shear 
viscosity. 

Resonant absorption and phase mixing are two inteiTe- 
lated mechanisms that may be involved in the energy cas¬ 
cade to small scales. Resonant absorption is an ideal pro¬ 
cess due to naturally occurring plasma and/or magnetic field 
inhomogeneities. In solar atmospheric plasmas, heating in¬ 
volving th e process of resonant absorption was first sug¬ 
gested by llonsoni (Il978h and has been extensively investi¬ 
gated in the literature afterwards. The theory of resonant 
MHD waves in the so lar atmosphere has been reviewed by 
IGoossens et alJ (1201 ll) . Resonant absorption has a strong im¬ 
pact on the temporal evolution of MHD kink waves in nonuni¬ 
form flux tubes because it produces a radial flux of wave en¬ 
ergy towards the nonuniform boundary of the flux tube (see 


e.g., iTataronisI II 975 I: l Abels-van Maanen & Weenink 
Poedts et al. 1989l:~ ^^egui et al.n201 iriPascoe et al 


1979; 


2013 


Goossens et alJ 1201 3h . The net result of the radial flux of 


energy is the simultaneous damping of the global kink mo¬ 
tion and the growing of rotational motions in the nonuni¬ 
form boundary of the tube. Although there is no direct ob¬ 
servational evidence for this mechanism in the solar plasma, 
the process has been checked by time-dependent numeri¬ 
cal simulations in which the MHD equations are evolved in 


time (see, e.g.. 

Terradas et alJl2006ll2008tlPascoe et alJl2013t 

IGoossens et al. 

120141). The numerical simulations show how 












































































2 


the energy of the global kink motion is transferred to the 
nonuniform boundary of the tube where small spatial scales 
are generated by phase mixing ( s ee. e.g..lPritchett & DawsonI 
[T 97 I [H^aerts & Priest] [1981 IcSblfloMh . The process 
of phase mixing occurs by the fact that the plasma is trans¬ 
versely nonuniform to the magnetic field direction, and so 
the Alfven frequency is spatially dependent. Alfven waves 
propagating in adjacent magnetic surfaces become more and 
more out of phase as time progresses. The consequence 
of this mixing of phases is the continuous decrease of the 
wave perturbations spatial scale across the m agnetic field 
(e.g.. lHevvaerts & Prieslill98^lMann et alJl995l) . In this sce¬ 
nario, the process of phase mixing would hypothetically keep 
working until the generated spatial scales become sufficiently 
small for Ohmic diffusion and/or viscosity to be operative. 
Then, efficient plasma heating due to Ohmic/v iscous dissi¬ 
pation of wave energy can take place (see, e.g.. lPoedts et alJ 
[ 1 ^[I^[T^ lOfman et al.|[l9^ . 

The goal of this paper is to investigate the time-dependent 
behaviour of MHD kink waves in nonuniform solar flux tubes 
by focussing on the processes of phase mixing and build¬ 
ing up of small scales. Therefore, we do not deal with the 
process of plasma heating itself, but we are interested in the 
previous step where the energy cascade form large scales to 
small scales takes place. Here we use a modal expansion 
to analytically express the MHD kink wave as a superposi¬ 
tion of discretized modes of the Alfven continuum. The full 
time-dependent evolution of MHD kink waves in nonuniform 
tubes has been investigated in the past by solving the initial- 
value problem analytically using the Laplace tra nsform (see 
iLee & Robertsl[T9^[RudWman & Robertsll2002^ and from a 

purely numerical point of view (see. e.g.. lTerradas et alj2006l 

[BollPascoe'eralMl among others). However, in the 
present paper we choose to follow a semi-analytic process 
that has some advantages compared to full numerical sim¬ 
ulations and is much simpler than the analytic method in¬ 
volving the Laplace transform. On the one hand, the ana¬ 
lytic part of the method allows an in-depth understanding of 
the physics behind the resonant absorption and phase mix¬ 
ing processes. On the other hand, since the numerical part 
of the method essentially consists in evaluating analytic ex¬ 
pressions, it is much faster than full numerical simulations 
and is free of the inherent artificial dissipation and resolution 
limitation of the numerical codes. The theory used in t his pa¬ 
per is inspired and follows closely that developed by ICallvl 
(119911) . In essence, we perfor m an extension to cylindrical 
geometry of the planar case of ICallvl (1199 Ih . The method 
has also been used in a number of other previous papers 


(see, e.g.. ICallv & Sedlaceklll992l 

1994 iMann etalJll99.5t 

ICallv & Maddisonlll997HTrrrv et al. 

Il997^ but. to the best of 


our knowledge, this method has never been used before to 
study the temporal evolution of MHD kink waves in nonuni¬ 
form cylindrical flux tubes. 

This paper is organized as follows. Sectioncontains the 
equilibrium model and the basic equations. The mathemati¬ 
cal method is explained in Section!?] The process of resonant 
absorption and the generation of small spatial scales due to 
phase mixing are investigated in Section|4| Later, Section|5]is 
devoted to the damping of the global kink motion and a com¬ 
parison to the usual global quasi-mode approach is performed. 
Finally, a discussion on the results and their implications are 
given in Section|6| 



Fig. 1 .— Slcetch of the cylindrical flux tube model used in this work. 

2. MODEL AND BASIC EQUATIONS 
2.1. Equilibrium 

We consider a simple model for a magnetic flux tube in the 
solar atmosphere. The sketch of the model is displayed in Fig- 
ure[T] The equilibrium configuration is composed of a straight 
magnetic cylinder of radius R embedded in a uniform and in¬ 
finite plasma. A cylindrical coordinate system is used, with 
r, ip, and z representing the radial, azimuthal, and longitudi¬ 
nal coordinates, respectively. The magnetic held is straight 
and along the axis of the cylinder, namely B = Blj, with B 
constant. The thermal pressure is constant everywhere. The 
density, p, is uniform in the azimuthal and longitudinal direc¬ 
tions, but nonuniform in the radial direction, namely 

f Pi, if r < n, 

Pir) = j Ptr(T), if ri < r < r 2 , (2) 

( Pe, if r > r2, 

with 

ri=R-^, r2 = R+^. (3) 

We use the subscripts ‘i’ and ‘e’ to denote quantities related 
to the internal and external plasmas, respectively, and the sub¬ 
script ‘tr’ to denote quantities related to the transitional layer. 
In Equation (|2|i, pi and pe are the internal and external con¬ 
stant densities. We set pi > pe to represent an overdense 
tube. In turn, ptr(T) denotes a nonuniform density profile 
that continuously connects the internal plasma to the exter¬ 
nal plasma across a layer of thickness / 6 [0,2/?]. The limit 
/ = 0 corresponds to a tube with a piecewise constant density 
(ri - r 2 - R), while the limit / = 2R corresponds to a fully 
nonuniform tube in the radial direction (ri = 0 and r 2 = 2R). 

2.2. Linear incompressible MHD waves 

We investigate linear MHD waves superimposed on the 
equilibrium ma gnetic flux tube. We use the idealMHD equa¬ 
tions (see, e.g.. lGoedbloed & Poedtsll2004tlPriesdl2014l) . We 
assume small perturbations over the static background plasma 
and linearize the ideal MHD equations, i.e., we only keep up 
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to linear terms in the perturbations. The resulting governing to Using Equation (| 6 ]l and combining Equations (fTTT i and 
equations are dll, we find that the redefined and are related to as 


^W§=-Vp' + -(VxB')xB, (4) 

B' = Vx(^xB), (5) 

where ^ is the plasma Lagrangian displacement, B' is the 
magnetic field Eulerian perturbation, p' is the thermal or gas 
pressure Eulerian perturbation, and p is the magnetic permit¬ 
tivity. Eor simplicity, we shall consider in this work that 
the plasma perturbations are incompressible. Hence, Equa¬ 
tions (IHl and (l5]l are supplemented with the incompressibility 
condition, namely 

V-^ = 0 . ( 6 ) 

We use the incompressibility condition and some vector iden¬ 
tities to rewrite Equation (Q as 

B' = B ■ V^. (7) 


■ 


nijr ld{r^r) 


kj -I- jr^ r dr 
kz 1 d (r^r) 


kj + rrP- jr'^ r dr 


In turn, the equation that relates P' with is 


P' = 


1 


P YfP jyl 


d^ _ B\2\ld(r^ 

dr 




(13) 

(14) 

(15) 


We combine Equations (fTOl i and (fTSl) to eliminate P' and ob¬ 
tain an equation involving only, namely 


Equation (l7]i allows us to elimate the magnetic field pertur¬ 
bation from the equations. We define the total (thermal 4 
magnetic) pressure perturbation as P' = p' 4 - B ■ B 


compact form, namely 

= -VP', 

where X.a is the Alfven operator given by 

52 g2 g2 


where the Alfven wave operator, X.A, and the surface wave 
operator, Xs, are 


in a 

d 2 

( 8 ) 


m 2 \ 52 

P. j Qy2 

(9) 

- 



ki + 


rdr 


+ ki- 


(17) 


(18) 


Eor later use, we explicitly write the three components of 
Equation ([ 8 ]l, namely 


d^ 


p{r) 


dfi 


B2 52 ^ 

1 dP’ 

p dz^] 


B2 52 \ 

1 dP' 

|<N 

1 ^ 

1 ^ 

r 0 (p 

B2 52 ^ 

i dP' 

p dz^i 

II 

1 


( 10 ) 

( 11 ) 

( 12 ) 


where and denote the radial, azimuthal, and longi¬ 

tudinal components of the Lagrangian displacement, respec¬ 
tively. 

Since the equilibrium is uniform in the ip- and z-directions, 
we can restrict ourselves to study the individual Eourier com¬ 
ponents of the perturbations along these directions. Different 
Eourier components do not interact with each other in the lin¬ 
ear regime. Hence the perturbations are put proportional to 
exp(/m^ 4 ikzZ), where m and kz are the azimuthal and longi¬ 
tudinal wavenumbers, respectively. Only integer values of m 
are possible. Kink waves are characterized by m = 1. These 
are the only waves that can displace the axis of the flux tube 
and can move it as a whole. Concerning the z-dependence, 
note that waves can be either standing or propagating, but the 
two types are equivalent from the mathematical point of view. 
Erom here on, we only retain the explicit dependence of the 
perturbations on the radial coordinate and time. 

We use as our main variable. Eor convenience, we rede¬ 
fine and ^z so that they incorporate the imaginary number 
i, namely and i^z We do so to avoid the pres¬ 

ence of imaginary numbers in the equations. Physically, the 
factor i accounts for a phase difference of 7r/2 with respect 


Equation (fTSl l governs incompressible MHD waves on the 
nonuniform flux tube and is the main equation of this investi¬ 
gation. 

When the density is uniform, dp{r)ldr - 0 and the second 
term on the left-hand side of Equation ( fThl l vanishes. The re¬ 
sulting equation is X.a-Cs& - 0. This equation has two types 
of solutions, namely classic Alfven waves and surface MHD 
waves, which are decoupled when the density is uniform and 
so they can be studie d separately. Readers are referred to 
iGoossens et al.l (120121) where the properties of both solutions 
were extensively discussed. On the one hand, the solutions as¬ 
sociated with the Alfven wave operator, i.e., the solutions to 
X.A& - 0 cortespond to classic Alfven waves. Classic Alfven 
modes have no motions along the magnetic field direction and 
have P' = 0. On the other hand, the solutions associated with 
the surface wave operator, i.e., the solutions to X.s^r - 0 cor¬ 
respond to surface MHD waves that have P’ + 0. Eor those 
modes, P' is maximum at the boundary of the tube and decays 
far away from the tube. 

Conversely, when the density is nonuniform dp{r')jdr + 0 
and the full Equation (fTSl ) must be considered. In such a 
case, classic Alfven waves and surface MHD waves are no 
longer separate solutions. We cannot clearly distinguish be¬ 
tween classic Alfven waves and surfa ce MHD waves (see a n 
extensive discussion on this matter in IGoossens et al.l| 2012 h . 
The two types of solutions get coupled due to plasma inho¬ 
mogeneity, whose effect is present in the second term on the 
left-hand side of Equation ( fThl l. Therefore, the solutions to the 
full Equation ( fTSl ) have mixed properties in the region with 
nonuniform density, i.e., in the nonuniform boundary layer of 
the flux tube. In this work we shall see that this fact has a 
strong impact on the temporal evol ution of the plasma mo¬ 
tions (see also IGoossens et ani2014l) . 
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As stated before, we assume that the MHD waves are in¬ 
compressible. The condition of incompressibility is used be¬ 
cause it greatly simplifies the mathematical analysis and is 
adequate for the type of transverse MHD waves studied in 
this paper. It has been shown that in the long-wavelength 
limit the fundamental radial mode of the MHD kink wave 
is almost incompressible to a high degree of accuracy (see 
iGoossens et al.ll2Ci0^l2012|) . Compressibility is of OiRIXf', 
where A = 2jilk^ is the longitudinal wavelength. The behav¬ 
ior of the fundamental radial mode is essentially the same in 
both compressible and incompressible cases when RfA <si 1 
or, equivalently, k^R <& 1. The condition of incompressibil¬ 
ity also imposes several limitations to our study that must be 
fairly acknowledged. Due to the condition of incompress- 
ibility there are no bo dy MHD waves in the system (see 
lEdwin & Robertslll983L for the distinction between surface 
and body waves). Body modes produce compression and 
rarefaction of the plasma in the flux tube and are not in¬ 
cluded in the analysis. However, this i s not a very strong 
restriction for our work. IGoossens et al.l (1201 2 |) showed that 
the fundamental radial mode of the MHD kink wave has the 
typical properties of a surface wave in the long-wavelength 
limit. Als o, radial harmonics a re absent in the incompress¬ 
ible limit (IGoossens et alJlT012l) . and lea ky waves in the ex¬ 
ternal plasma are not possible either (see ICallMll986h . Both 
radial harmonics and leaky waves are important in the tran¬ 
sitory pha se that follows the e xcitation of standing oscilla¬ 
tions (see iTerradas et al1l2006ll and are also involved in the 
propagation of sh ort-wavelength waves along the tube (see 
lOliver et^l2014t) . Hence, to be consistent with the incom¬ 
pressibility assumption, we shall restrict ourselves to study 
the fundamental radial mode of long-wavelength MHD kink 
waves. 


internal plasma is 

^r{r, t) = Ai(f)C , if r < ri, (20) 

where Ai(f) is a time-dependent amplitude. In turn, in the 
external plasma, i.e., for r > r 2 , we require to vanish when 
r ^ oo, so that solution in the external plasma is 

^r(r, t) = As(t)K'^ ik^r ), if r > ra, ( 21 ) 

where again Ae(f) is a time-dependent amplitude. 

We recall that we restrict our analysis to long-wavelength 
kink waves, i.e., we assume k^R « 1 and set m = 1. To 
simplify matters, we perform asymptotic expansions of and 
K'm for k,R <si 1 , i.e., for small arguments, and keep the first 
term in the expansions only. Therefore, for long-wavelength 
kink waves the form of in the internal and external plasmas 
simplifies to 

if r<ru ( 22 ) 

f,(r,f)^Ae(f)r-^ if r>r2, (23) 

where, with no loss of generality, all constants have been ab¬ 
sorbed into Ai(f) and Ae(f)- 

3.2. Solution in the nonuniform boundary 

In the nonuniform layer, i.e., for ri < r < r 2 , there is a 
continuum of Alfven modes that couple to the surface MHD 
kink wave. We cannot remove Alfven waves from the anal¬ 
ysis, and this fact has important physical consequences. We 
must necess arily consider the full Equation (fTbl l. 

Eollowing lCa^llvl (119911) . we preform a modal expansion of 
the displacement in the nonuniform boundary of the flux tube. 
We write ^^(r, f) as a generalized Eourier series, namely 


3. MATHEMATICAL METHOD 

Here we find the solutions to Equation (fTbl l. We stress 
that we do not assume that there exists a global mode in the 
flux tube with a certain frequency. Instead, we retain the full 
temporal dependence of the perturbations. The solutions to 
Equation ( fTSb in the regions of the flux tube with uniform and 
nonuniform density are analyzed separately. 


3.1. Solution in the uniform internal and external regions 

As discussed before, classic Alfven waves and surface 
MHD waves are decoupled when the density in uniform. We 
are interested in studying the MHD kink wave, which is a 
surface mode in the incompressible limit. So, for the sake of 
simplicity we remove the classic Alfven waves from the scene 
in the regions with uniform density. Then, Equation (fTbl) can 
be cast as 




dr^ 


K + 


3m^ 


/ 2 

ill 

1 

\k^ + —\ + 

k'^ -^ 

— 

[r 


yl 


1 ^ 
r dr 


(19) 


Equation (fT^ is a differ ential equation related to the modi- 
hed Bessel equation (see lAbramowitz & Stegunlll97^ . The 
solution to Equation ( fT9l l is a linear combination of functions 
I'„fkj) and Kl„{k^r), where 4 b and are the usual modihed 
Bessel functions and the prime ' denotes a derivative of the 
modihed Bessel function with respect to its argument. 

In the internal plasma, i.e., for r < ri, we require fr to be 
regular at r = 0. Hence, the solution to Equation ( fT9] l in the 


&ir, 0 = ^ an(t)if/„(r), (24) 

n=l 


where a„{t) is the time-dependent amplitude of the «-th base 
function i4n(r). As appropriate in a cylindrical coordinate sys¬ 
tem, we choose functions i 4 n(r) to be eigenfunctions of the 
regular Sturm-Liouville system dehned by the Bessel differ¬ 
ential equation 


dftf/ 1 dtf 

d^2 f. 



(25) 


together with the boundary conditions 


dlA y. 

_= 0 , at r = ri, 
dr 

(26) 

2 dd/ 

—if H-= 0 , at r - r 2 . 

r dr 

(27) 


These boundary conditions are chosen for the eigenfunctions 
to match the spatial behavior of at r = ri and r - r 2 (see 
Section ITTI) . The eigenfunctions, have the following 

general form, 

^n(r) = Qn [Jx{Anr)Y[{Anr,) - J[{Anr,)Yf Anr)] , (28) 


where A„ is the «-th eigenvalue, Ji and Yi are the usual Bessel 
functions of order 1 , and Qn is a constant determined by 
imposing the orthogonality condition on the eigenfunctions, 
namely 

1 n 

- J i/r„(r)i/r„fr)rdr = (29) 
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In turn, the eigenvalues, A„, are the solutions of the transcen¬ 
dental equation, 

Ar2 Ji(Ar2)Y;(An) - j;(An)Yi(Ar2) 

2 ~ j;(Ar2)Y;(An) - j;(An)Y;(Ar2y 

The constants Q„ and the eigenvalues An can be numerically 
obtained. 

Figure |2] displays the eigenfunctions of the first four gen¬ 
eralized Fourier modes. The larger the order of the Fourier 
mode, the more spatial oscillations of the eigenfunction. In 
this formalism, large spatial scales are related to Fourier 
modes of low order, while small spatial scales are related to 
Fourier modes of high order. Therefore, Fourier modes of 
low/high order are expected to lose/gain weight as time pro¬ 
gresses during t he processes of phase mixing and build ing up 
of small scales (lCallvllT9^ iCallv & Maddisonlll997h . We 
shall confirm this in Section l4~^ 

Next, we substitute Equation (l24l i into Equation (fThl) and 
arrive at a relation for the time-dependent amplitudes, a„(t), 
namely 


p(r)-L&ilfnir) + [k] 


w^\/ di/rn(r) 
j\ dr 


+kl — £s^n(r)a„(t) = 0 , 

P 



d^a «(0 


(31) 


for n = 1,2,.... Next, we set the temporal dependence of 
as exp(-ia>t), multiply Equation OTI) by i/'n'(r), and integrate 
the resulting equation over the interval [ri, r 2 ] to eliminate the 
dependence on r. Then, a generalized eigenvalue problem can 
be cast, namely 

Ha^w^Ma, (32) 

where a = [ai,a 2 , ■ ■ -V is the right eigenvector, ci/ is the 
eigenvalue, and H and M are square matrices whose elements 
are 




^ p I X 


t/r„(r)£si/r„,(r)rdr. 


1 

iWn«' = yJ [pir)£s^/fn'{r) 


dp(r) 


, m 


dipn'ir) , 1 , , , 

-IT 


(33) 


iA„(r)rdr. 

(34) 


When the density is uniform, all the eigenvalues are degen¬ 
erate, namely = a>\, where a>A - k^VA is the Alfven 
frequency, with va - Bj yfjlp the Alfven velocity. When 
the density is nonuniform the Alfven frequency is spatially- 
dependent, i.e., a continuum of Alfven frequencies is present. 
In such a case, it can be shown that the eigenvalues are neces¬ 
sarily in the interval ^ g, where UA,i - and 

WA.e = KvAfi are the intern al and extern al Alfven frequencies, 
respectively (see details in lCallviri991h . 

Once the eigenvalue problem is solved, the time-dependent 
amplitudes can be expressed as a superposition of eigen- 
modes, namely 


CO 

««(0 - {c„' COS (w„-f) + dn' sin (w„-f)], (35) 

w' = l 


where is the n'-th eigenvalue, is the n-th component 



Fig. 2.— Eigenfunctions of first four generalized Fourier modes in a nonuni¬ 
form layer with IjR = 1. For representation purposes, the eigenfunctions ai‘e 
normalized so that ^„(ri) = 1. 


of the n'-th right eigenvector, and c = [ci,C 2 ,.. .]^ and d = 
{d\,d 2 , ■ ■ depend upon the initial conditions as 


c=y 8 -'a(f = 0 ), 


d=y 8 * 



(36) 

(37) 


where j8 is a matrix whose elements are y6„„/. The components 
of a and their derivatives at f = 0 are computed from the initial 
conditions for as 


Unit 

dan 

dr 


-0>-f 

Jri 

-r 

t=0 Uri 


^r(r, t = Q)i//„(r)rdr, 


d^r(r, t) 


dt 


{r,t=0) 


ipn(r)rdr. 


(38) 

(39) 


Finally, the full expression of t) in the nonuniform layer 
is 


CO oo 


nn' [Cn' COS {oJn't) + dn' sin (w„'f)] ^n{r). 

n=\ n'—\ 

(40) 


3.3. Alfven continuum modes 

The generalized eigenvalue problem of Equation (l32l i in¬ 
volve infinite matrices. To solve the eigenvalue problem nu¬ 
merically, we must truncate the generalized Eourier series so 
that only the first N Eourier modes are considered. Hence, we 
must replace oo by A in the upper limit of all summations. 
Then, the infinite matrices become finite N xN matrices and, 
consequently, the number of eigenvalues is N. Eurthermore, 
in order to shed light on the physics behind Equation (l40l) . we 
recast it in the following alternative form, 

N 

fr(r, 0 = ^ [cn COS (w„f) H- dn sin (w„f)] (f)„(r), (41) 


with 


N 

( 0 - 

«' = 1 


(42) 
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Fig. 3.— Left: Alfven contiuum of frequencies (dotted line) as function of position in the nonuniform layer. The superimposed crosses correspond to the discrete 
eigenfrequencies for N = 31. Right: Eigenfunctions of three selected discrete modes for N = 31 (solid lines) and for N = 201 (dotted lines). In the computations 
we used a sinusoidal transition of density, p jpe = 5, l/R = 1, m = 1, and k-R = ;r/100. 


As explained by ICallvl (1199 Ih . functions <f>„(r), and their as¬ 
sociated frequencies aj„, play the role of the Alfven contiuum 
modes. The whole set of these N modes represents a dis¬ 
cretized version of the Alfven contiuum. Equation (HTt is use¬ 
ful because it allows us to physically understand that the to¬ 
tal displacement, t), is the superposition of the displace¬ 
ments produced by all the modes of the Alfven continuum. 
These discrete ‘Alfve n contiuum mo des’ are equivalent to the 
modes discussed by iBarstonI (119641) in the context of elec- 
trostatic oscill ations in a nonuniform plasma. As shown by 
Barstonl (1 19641) . no global mode is present when the density is 
nonuniform. 

Figure [3 (left panel) displays the Alfven continuum of fre¬ 
quencies, ojpXr), as function of position in the nonuniform 
layer for a particular set of parameters given in the caption of 
the Figure. In this computation we used a sinusoidal variation 
of density in the nonuniform layer, namely 



We overplot in Figure |3] (left panel) the discrete eigenfreque- 
cies, w„, that are solutions to the truncated eigenvalue prob¬ 
lem using A = 31. We see that the Alfven contiuum is cor¬ 
rectly recovered by the discrete eigenfrequencies. We have 
selected three of these discrete modes and have plotted their 
eigenfunctions in the right panel of Figu re |3] This Figure can 
be compared to Figure 8 of ICallvl (Il991h for equivalent eigen¬ 
functions in planar geometry. The eigenfunctions are discon¬ 
tinuous (singular) at the specific position in the nonuniform 
layer where the discrete eigenfrequency matches the local 
Alfven frequency. We notice both logarithmic and step dis¬ 
continuities in the plotted eigenfunctions, as consistent with 
the expec ted behavior of the Alfven continuum modes (see 
details in iGoedbloed & Poedtsll200^ . The singularities are 
not well resolved in the right panel of Figure [3] because, on 
purpose, we considered a relatively small value of N. As a 
consequence, there are too few regular Fourier modes to ac¬ 
curately capture the singular behavior of the eigenfunctions. 
Also note the presence of wobbles as a sort of Gibbs phe¬ 
nomenon in the eigenfunctions. These issues can easily be 
overcome by increasing the value of N. The larger the value 
of A, the larger the number of Fourier modes and so the bet¬ 


ter representation of the contiuum modes. For instance, we 
have increased the number of Fourier modes to A = 201 
and have computed again the eigenfunctions of the discrete 
‘Alfven contiuum modes’ whose frequencies are the same as 
those displayed in Figure [3] (right panel). We have overplot¬ 
ted these new eigenfunctions using dotted lines. The singular 
behavior of the eigenfunctions is now well captured and the 
wobbles are much less noticeable than for A = 31. 

Note in Figure [3 (right panel) that the ‘wings’ of the 
eigenfunctions penetrate in the internal and external plasmas. 
Hence, in the regions with uniform density the added contri¬ 
bution of the ‘wings’ of all the continuum modes are respon¬ 
sible for the total displacement. In addition, the fact that ev¬ 
ery continuum mode has a different frequency has the conse¬ 
quence that the oscillations in the nonuniform layer will even¬ 
tually get out of phase as time progress es. In essence, this 
is the process of phase mixing (see, e.g.. lPritchett & DawsonI 
ll978t[Hevvaerts & Prieslill983l) . A question then arises: what 
happens to the displacement in the internal and external re¬ 
gions with uniform density as phase mixing occurs in the 
nonuniform layer? We anticipate that the net result is the ap¬ 
parent damping of the global kink motion. This is explored in 
detail in SectionH] 


4. TEMPORAL EVOLUTION OF THE DISPLACEMENT 

Here, we show the result of the temporal evolution of the 
displacement components in the flux tube using the modal ex¬ 
pansion method. As initial condition for the radial component 
of the displacement we set 


^Ar, f = 0) = < 


f^o, 

^0 


e ( r2'A 


if r < r\, 
if ri < r < r 2 , 

if r > r 2 . 


(44) 


and t= 0 ) - 0 for r, where is the arbitrary internal 

amplitude at f = 0. The azimuthal, and longitudinal, 
components of the displacement depend on through Equa¬ 
tions ( fT3] ) and (fT4l i. respectively. With this initial condition, 
only the first generalized Fourier mode, i.e., the mode with 
the largest spatial scale, contributes to the total displacement 
at f = 0. We choose this initial condition on purpose to investi¬ 
gate the process of building up of small scales that are related 
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Fig. 4. — Normalized amplitudes of the discrete ‘Alfven contiuum modes’ as 
func tion of their frequencies according to the initial condition given in Equa¬ 
tion I44> . The red diamonds correspond to IjR = 0.2 and the blue crosses are 
for IjR = 1. The vertical dashed lines denote the corresponding quasi-mode 
frequencies. 

to Fourier modes of large order. In all computations we use 

= 201 to assure a good accuracy of the modal scheme. Un¬ 
less otherwise stated, we consider a sinusoidal transition of 
density (Equation (l43T l) and use the following set of parame¬ 
ters: pi/pe - 5,m - 1, and k^R = ;7r/100. 

First of all, we solve the generalized eigenvalue problem 
(Equation (13^ 1 for two values of the nonuniform layer thick¬ 
ness, namely a thin layer with l/R = 0.2 and a thick layer with 
l/R - 1. Then, we compute the amplitudes, c„, of the discrete 
‘Alfven contiuum modes’ using the initial condition given in 
Equation (l44l) . Note that d„ - 0 since 
The normalized values of c„ are plotted against their corre¬ 
sponding eigenfrequencies in Figure 0] In the case of a thin 
layer, we find a sharp peak around a specific frequency, which 
indicates that the ‘Alfven contiuum modes’ that mostly con¬ 
tribute to the total displacement are those whose frequencies 
are near that specific frequency. However, in the case of a 
thick layer a broader spectrum of modes contributes to the to¬ 
tal displacement. We shall discuss in Section |5] the relation 
between the position and width of the peak in the frequency 
spectrum and the so-called quasi-mode freq uency and damp¬ 
ing rate (see also lAndries & Goossensl[2007n . 

Next, the temporal evolution of the displacement compo¬ 
nents is computed. Figures |5] and |6] display the temporal 
evolution for l/R — 0.2 (Figure |5]l and l/R - 1 (Figure |6]l. 
See also the accompanying animation^ Time is expressed in 
units of Pk - iTi/cJk, where Pk is the kink mode period of a 
homoge neous thin tube, with u)k the so-called kink frequency 
(see, e.g. lEdwin & Robertsll 19831) given by 


Oik 


Pi^L +P="i,e 


Pi +Pe 


(45) 


In the following subsections, we discuss in detail the dynam¬ 
ics shown in Figures |5] and |6] These results can be compared 
to those obtaine d from full numerical simu lations of kink os¬ 
cillations (e.g.. iTerradas et akl 120061 l2008h . Also, although 
found in a very different context, our results show many sim¬ 


* Movies can be downloaded from 

http://www.uib.es/depart/dfs/Solar/movies_SolerTerradas2815.zip 


ilarities to those obtained bv iMann et al.l (119951) in the case of 
resonant waves in the magnetosphere. 

4.1. Resonant absorption and flux of energy to the 
nonuniform boundary 

The temporal evolution of the displacement components 
displayed in Figures |5] and |6] shows that, as time evolves, 
the amplitude of the oscillations in the nonuniform layer in¬ 
creases, while the amplitude in the internal and external plas¬ 
mas decreases. The global transverse oscillation damps and, 
after several periods, the oscillations are essentially polarized 
in the azimuthal direction. In other words, the initially trans¬ 
verse motion becomes rotational. As time increases, the az¬ 
imuthal component of the displacement in the nonuniform 
layer clearly dominates the dynamics, while both radial and 
longitudinal components are much smaller than the azimuthal 
component. This behavior is fully consistent with the pres¬ 
ence of an energy flux towards the nonunifo rm layer due to 
the pr ocess of resonant absorption (see, e.g., iGoossens et al.l 

IMl- 

To further illustrate the energy flux towards the nonuniform 
boundary of the tube due to resonant absorption, we compute 
the total (kin etic + magnetic) energy density, E, as (see, e.g., 
IWalkedl200l 

£= ijpjv'f+ (46) 

where v' = d^/dt and B' = Bd^/dz are the velocity and mag¬ 
netic field perturbations, respectively. Figure |7] (left panels) 
shows the evolution with time of the spatial distribution of the 
energy density for the oscillations displayed in Figures |5] and 
|6l As time increases, more and more energy is fed into the 
nonuniform layer. Eventually, the shape of the energy den¬ 
sity becomes a Gaussian-like function localized around the 
nonuniform boundary of the tube. An equiv alent result can 
be seen in Figure 8 of ITerradas et al.l (I2006^ ob tained from 
numer ical simulations, and also in Figure 7 of iMann et alJ 
(119951) for the case of magnetospheric waves. 

The right panels of Figure [T] display the integrated en¬ 
ergy as function of time in each region of the flux tube, 
namely the internal plasma, the nonuniform layer, and the 
external plasma. The amount of energy localized in the 
internal and external plasmas becomes negligible at suffi¬ 
ciently large times. The efficiency of resonant absorption 
increases as the nonuniform layer gets th i cker (see, e.g., 
Rae & RobertsI 119811: iLee & Roberts! 1 19861 : IGoossens et alJ 
19921: iRuderman & Robertsll2002^ . For this reason, the pro¬ 
cess of energy transfer is faster in the case of a thick layer 
than in the case of a thin layer. The energy flux towards the 
boundary of the tube is very fast when the tube is largely in¬ 
homogeneous. For instance, for l/R = 1 practically all the 
available energy is already localized around the boundary of 
the tube after two periods of the kink oscillation. The process 
is slower for l/R = 0.2, but the final result is the same: all the 
wave energy is transferred to the nonuniform boundary of the 
tube. 

4.2. Phase mixing and energy cascade to small scales 

Small spatial scales are generated in the inhomogeneous 
layer as the process of resonant absorption feeds energy into 
the nonuniform boundary of the tube. The temporal evolu¬ 
tion of the displacement clearly shows how these small spatial 
scales evolve due to phase mixing. Smaller and smaller scales 
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Fig. 5.— Temporal evolution of (left), (center), and (right) in a nonuniform tube with IjR = 0.2. Top, mid, and bottom panels are for tlP^ = 0, 3, and 
10, respectively. The vertical dashed lines denote the boundaries of the nonuniform layer. Note that the vertical scale is different in each panel. Movie can be 
downloaded from http://www.uib.es/depart/dfs/Solar/movies_SolerTerradas2015. zip 


are generated as time increases (s ee. e.g..lPritchett & DawsonI 
Il978t iHevvaerts & Prieslll983h . iMann et alJ (Il995h derived 
an approximate expression for the length scale generated due 
to phase mixing, Lp^, at a given position, r - ro, and as a 
function of time. The expression is 

“ 77 77“i 7’ 

\da)Aldrl^ro t 

where \d(jjfi,ldr\r^r^ is the absolute value of the gradient of the 
Alfven frequency at r = tq. Figure |8] (left panel) displays Z™ 
at r = /? as a function of time computed from Equation (l47t . 
In that graph, we overplot using symbols the length scale di¬ 
rectly estimated from the evolution of the displacement ob¬ 
tained from the modal expansion method. To do that, we have 
taken the azimuthal component of the displacement and, for a 
fixed time, have measured the distance between the two con¬ 
secutive extrema closer to r - R. We find a remarkable agree¬ 
ment between the length scale estimated by this method and 
Equation ( |47] ). As expected, the phase-mixing length scale 
gets shorter as time increases. We also see that the thinner the 
nonuniform layer, the shorter Lp^. This last result is consis¬ 
tent with the fact that the density profile is steeper for a thin 
layer than for a thick layer. Hence, the thinner the nonuniform 


layer, the larger the gradient of the Alfven frequency, and so 
the smaller the right-hand side of Equation (l47l i for a fixed 
time. 

In the modal expansion, different spatial scales are associ¬ 
ated to Eourier modes of different order. Long/short scales are 
associated to Eourier modes of low/high order. By substitut¬ 
ing the Eourier expansion of the displacement (Equation (l24l ll 
into the expression for the energy density (Equation (l46l) I. it 
can be shown that the contribution of the n-th Eourier mode 
to the total energy is proportional to a^. Eigure|8] (right panel) 
displays the temporal evolution of for n = 1, 5, 10, and 20 
in a flux tube with l/R - 1. This Eigure evidences the energy 
cascade from large scales to small scales due to phase mixing. 
Eourier modes of high order (i.e., small scales) gain weight in 
the total energy as time increases, while Eourier modes of low 
order (i.e., large scales) lose weight. As time evolves, energy 
is distributed among smaller and smaller scales. 

The process of phase mixing works indefinitely in our 
model since there is no dissipative mechanism included. On 
purpose we used the ideal MHD equations to observe the pro¬ 
cess of energy cascade. Let us consider the scenario in which 
magnetic diffusion is included. We can relate Lpm with the 
typical length scale, Iq, that appears in the definition of the 
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Fig. 6. — Same as Figure|^but with//i? = 1. Movie can be downloaded from http://www.uib. es/depart/dfs/Solar/movies_SolerTerradas2815. zip 


magnetic Reynolds number (Equation ([T])). Then, it can be 
seen that /?,„ ~ owing to the effect of phase mixing, i.e., 
magnetic diffusion becomes more and more efhcient as time 
increases. In reality, the behavior of the MHD kink waves 
would be essentially the same as in the ideal case studied 
here until the generated spatial scales due to phase mixing 
are sufhciently small for d iffusion to become efhcient (see 
iGoedbloed & Poedtsll2004t) . Then, wave energy can be de¬ 
posited in the background medium and so heat the plasma 
via Ohmic/viscous heatin g (see, e.g.. lPoedts et al.ll989lll990l 
Il994t lOfman et al]ll995ll . An equivalent discussion is given 
by Mann et alJ d 19951) for the case of dissipation in the mag¬ 
netosphere due to Pedersen conductivity. 

5. DAMPING OF THE GLOBAL KINK MOTION AND COMPARISON 
TO THE QUASI-MODE 

Here, we focus on the damping of the global transverse 
motion and analyze how the amplitude of the displacement 
changes with time near the axis of the flux tube. The pur¬ 
pose of this Section is to compare the results from the modal 
expansion method, in which there are no global modes, with 
those of the usual global quasi-mode approach. 

The usual procedure to theoretically investigate the damp¬ 
ing of MHD kink waves in a transversely nonuniform flux 
tube is based on the assumption that there exists a global 


kink mode. In other words, it is assumed that a nonuniform 
tube supports a global kink mode that is the direct descen¬ 
dant of_the_^in^jnodeofan equivalent uniform tube (see, 
e.g.. lEdwin & Robertslll98^ . This hypothetical global kink 
mode represents a coordinated transverse motion of the flux 
tube, so that the whole plasma is as sumed to oscillate at the 
same frequency (for more details see lTirrv & Goossensll 19961 
and references therein). The frequency of the global kink 
mode is complex and the imaginary part of the frequency is 


Sedlacekll971l;lTataronis & Grossmannll973l: 

Rae & Roberts 

1981 

: Rael 1982HGoossens et al.l 19921: iRuderman & Roberts 

2002 

: ISoler et alJ 120131). However, complex eigenvalues 


are n o t possible in ideal MHD ( see, e.g., iPoedts & Kernel! 

i). 


119911 IGoedbloed & PoedtsI 120041) . A complex frequency 
cannot possibly correspond to a true ideal normal mode 
of the flux tube. Eor this reason, the complex mode that 
physically represents a damped global motion of the flux 
tube is called a ‘quasi-m ode’ or a ‘virtual mode’ (see also 
ISedlacek & Roberdl989h . 

The quasi-mode temporal dependence is proportional to 


cos (faiqmf) exp i-jt), where and y are the quasi-mode fre¬ 
quency and damping rate, respectively. Here we compute the 
quasi-mode frequency and damping rate with the method de- 































































































E/max(E) E/max(E) 


10 


l/R = 0.2 




l/R = 1 




Fig. 7.— Energy density as function of position at three different times (left) and integrated energy in each region of the flux tube as a function of time (right) 
in a nonuniform tube with IjR = 0.2 (top) and IjR = 1 (bottom). The remaining parameters are the same as in Figure[5] 




Fig. 8.— Left: Small spatial scale generated due to phase mixing at r « 7?. The solid line corresponds to Equation ED and the symbols are the spatial scale 
directly estimated from the evolution of the displacement in the modal expansion method. Right: Energy cascade from large scales to small scales. Temporal 
evolution of (in logarithmic scale) for n = I, 5, 10, and 20 in a flux tube with l/R = 1. 
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l/R = 0.2 l/R = 1 




Fig. 9.— Temporal evolution of ^r/^o at r = 0 in a nonuniform tube with IjR = 0.2 (left) and IjR = 1 (right). The blac k solid line is the result from the modal 
expansion method and the blue dashed line is the dependence predicted by the quasi-mode (computed with the method of ISoler et al.l2Ql3) . 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 


l/R = 1 



Fig. 10.— Power spectrum (normalized power vs. frequency) of at r = 0 displayed in Figure[9] The two vertical red dott ed lines denote th e internal and 
external Alfven frequencies, while the vertical blue dashed line denotes the quasi-mode frequency (computed with the method of ISoler et alJ201^ . 
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Fig. 11. — Temporal evolution of ^r/^o at r = 0 in a nonuniform tube with 
IjR = 1 and a very low density contrast of pi/pe = 1-2. The vertical dotted 
line denotes the switch fro m Gaussian-like to ex ponential-like envelopes ac¬ 
cording to Equation (10) of IPascoe et alj <20131) . The decayless phase takes 
place for tjPk > 20. 


scribed in ISoler et alJ ( l2013h . In short, we solve the singular 
MHD equations in the nonuniform layer using the Method of 
Frobeni us and find a ‘d ispersion relation’ for the quasi-mode 
(see also lHollweglll99Cil) . The numerical solution to the quasi¬ 
mode ‘dispersion relation’ provides us with Wqm and y. For 
the sake of simplicity we omit t he details of the q uasi-mode 
computation and refer readers to lSoler et alJ (1201 3|) for an ex¬ 
tensive explanatio n of the procedure . Importantly, we note 
that the method of ISoler et all (1201 3h is general and uses no 
approximation to compute the quasi-mode and, therefore, it is 
a pplicable to thick nonuniform layers. 

lAndries & GoossensI (l2007l) investigated the continuous 
spectrum of leaky MHD modes in a slab. These authors 
showed that the spectral measure associated with the con¬ 
tinuous leaky spectrum is not a monotonic function of the 
frequency, but peaks appear in the spectral measure around 
specific frequencies (see their Figure 2). According to 
lAndries & GoossensI (120071) . these specific frequencies corre¬ 
spond to the discrete dam ped leaky modes described in the lit¬ 
erature (e.g.. lCair^ll986[) . while the characteristic time scale 
associated with the decay of those continuum leaky modes 
is determined by the width of the peak s. We can perform a 
clear parallelism between the results of lAndries & GoossensI 


clear p a 

dlOOl i 


I for leaky modes and those obtained here for resonantly 
damped quasi-modes. In Figure |4] we see that the peak in 
the frequency distribution of the ‘Alfven continuum modes’ is 
wider for IjR - \ than for IjR - 0.2. This result is consistent 
with the fact that the damping rate of the quasi- mode increases 
as the nonuniform layer gets thicker (see, e.g., Lee & RobertsI 
' ISoler et S! 20131). and i 


Il986t iRuderman & RobertsI 120021 ISoler et ^ 


IS 


equiva lent to the behaviour discussed bv lAndries & GoossensI 
(l2007h for the decay of leaky modes. We have overplotted 
in Figure |4] the quasi-mode frequency. When IjR — 0.2 the 
quasi-mode frequency matches perfectly the sharp peak in 
the distribution of the ‘Alfven continuum modes’. However, 
when IjR - \ the quasi-mode frequency is smaller than the 
central frequency of the peak. When IjR - 1 it is unclear 
how the quasi-mode frequency is related to the distribution of 
the continuum frequencies. This result may suggest that the 
quasi-mode frequency is a good approximation to the effec¬ 


tive frequency of the global oscillation when the transitional 
layer is thin, but it may become less representative when thick 
layers are considered. 

In order to explore the accuracy of the quasi-mode com¬ 
pared to the actual temporal evolution, we display in Figure|9] 
the temporal evolution of at the axis of the flux tube. Again 
we consider two values of the nonuniform layer thickness, 
namely a thin layer with IjR — 0.2 and a thick layer with 
IjR - 1. These results correspond to the same computations 
displayed in Figures |5] and |6] but now we focus on the be¬ 
havior at r = 0. For comparison, we overplot in Figure |9] the 
temporal dependence of the quasi-mode. The modal expan¬ 
sion method and the quasi-mode essentially provide the same 
result when the nonuniform layer is thin (see the left panel of 
Figure 0. In that case the amplitude of near the flux tube 
axis decreases exponentially at the quasi-mode damping rate. 
However, as expected, the results are more different when the 
nonuniform layer is thick (see right panel of Figure |9]l. The 
result from the modal expansion method shows two distinct 
phases. In the first phase, approximately oscillates at the 
quasi-mode frequency, but its amplitude does not decrease as 
fast as the quasi-mode predicts. In the second phase, os¬ 
cillates at a frequency different from that of the quasi-mode 
and its amplitude remains roughly constant, i.e., the oscilla¬ 
tions in this second phase are almost decayless. In the right 
panel of Figure|9]the switch from the first phase to the second 
phase takes place at tjPk ~ 3. Hence, the quasi-mode offers 
a poor approximation to the actual transverse displacement of 
the flux tube axis in the case of thick nonuniform layers and 
large times. 

To shed more light on the similarities and differences be¬ 
tween the results from the modal expansion method and those 
from the quasi-mode, we display in Figure[T0]the power spec¬ 
tra of the signals of Figure |9] In the case of the result for 
a thin nonuniform layer (left panel), we find that the power 
spectrum has a pronounced peak around the quasi-mode fre¬ 
quency. This result confirms that the flux tube axis oscil¬ 
lates with an effective frequency that corresponds to the quasi¬ 
mode frequency. However, contrary to the case for a thin 
layer, we do not find such a clear, sharp peak in the power 
spectrum for a thick layer (right panel). Although with some 
differences, in both cases the power spectrum at r = 0 re¬ 
sembles the frequency distribution of the ‘Alfven continuum 
modes’ in the nonuniform layer (compare Figuresl4land[T0ll. 

In the power spectrum for thick layers (Figure [TO] right 
panel), there is significant power around the quasi-mode fre¬ 
quency (as for thin layers) but we also find significant power 
in those frequencies near the internal Alfven frequency. An 
explanation for the shape of the power spectrum is found in 
the presence of the two phases of the oscillation discussed 
above (see again the right panel of Figure |9|l. The scenario 
is as follows. In the first phase of the oscillation, the flux 
tube axis oscillates at an effective frequency that matches the 
quasi-mode frequency. The oscillation amplitude decreases 
rapidly and then it starts the second phase with a roughly 
constant amplitude. The oscillation frequency in the second 
phase turns out to be very close to the internal Alfven fre¬ 
quency. Hence, both the quasi-mode frequency and the in¬ 
ternal Alfven frequency significantly contribute to the power 
spectrum when the nonuniform layer is thick. Conversely, no 
trace of the internal Alfven frequency is found in the power 
spectrum when the nonuniform layer is thin since the quasi¬ 
mode dominates for the whole duration of the temporal evo¬ 
lution displayed in Figure |9] (left panel). 
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Finally, it is interesting to test whether the modal expan¬ 
sion method is able to recover the initial Gaussian-like enve¬ 
lope of the displ acement obtained in time-dependent numer¬ 
ical simulations ([P ascoe et al.ll2012l 120131 iHood et ^120131 
Ruderman & Terradasll20i:^ The quasi-mode damping pro¬ 

file is purely exponential, but previous works showed that a 
Gaussian profile shows up in early stages of the oscillation if 
the density contrast is low. We show in Figure[TT]the temporal 
evolution of fr at r = 0 for IjR - \ and pi/pe = 1-2. The ini¬ 
tial Gaussia n-shaped envelope i s evident. According to Equa¬ 
tion (10) of iPascoe et akl ( 120131) . the switch from Gaussian to 
exponential envelopes occurs sltjPk ~ 11 for the parameters 
considered in Figure [TT] and agrees well with the result dis¬ 
played in that Figure. In addition, we again see in Figure [TT] 
the presence of the second, decayless phase of the oscillations 
when tjPk > 20. 

6. DISCUSSION 

In this paper, we have studied the temporal evolution of 
MHD ki nk waves in nonuniform magnetic flux tubes. In¬ 
spired bv ICallvl (1199 II) . we have used a semi-analytic method 
based on expressing the Lagrangian displacement of the in¬ 
compressible MHD kink wave as a superposition of dis¬ 
cretized modes of the Alfven continuum. The behavior of the 
MHD kink wave obtained with this method matches that seen 


in full numerical simulations (e.g.. 

Terradas et al.ll2006l 120081 

IPascoe et al.ll2013l;lGoossens et al. 

120141). Hence, the present 


paper shows an alternative method that can be used in future 
works to include the full dynamics of MHD kink waves in, 
e.g., seismology applications or forward modelling. 


The evolution of the MHD kink wave with time shows how 
the processes of resonant absorption and phase mixing oper¬ 
ate. Due to the combined effect of both mechanisms, the ini¬ 
tial global kink oscillation of the flux tube is eventually trans¬ 
formed into small-scale rotational motions in the nonuniform 
boundary of the tube. Resonant absorption transfers wave en¬ 
ergy towards the nonuniform boundary, where phase mixing 
is responsible for the energy cascade from large spatial scales 
to small spatial scales. The two processes occur simultane¬ 
ously and, indeed, they are intimately linked; both processes 
are caused by plasma and/or magnetic field inhomogeneity 
across the flux tube. Such inhomogeneities naturally occur in 
magnetic flux tubes of the solar atmosphere. 

Due to the use of the linear ideal MHD equations, our 
analysis misses the very last step of the evolution of the 
MHD kink wave when the spatial scales generated by phase 
mixing are sufficiently small: heating of the background 
plasma due to Ohmic/viscous dissipation of wave energy. To 
study the heating, the nonlinear dissipative MHD equations 
should be solv ed using, necessari l y, ful l numerical simula¬ 
tions (see, e.g.. iPoedts et al.l[T989L 119901 Il994t lOfman et ^ 
Il995h . In this work we have studied the evolution of the 
MHD kink wave as an initial-value problem. Hence, the to¬ 
tal energy available in the system is determined by the ini¬ 
tial condition and, therefore, the maximum amount of en¬ 
ergy that can be deposited in the plasma in the form of 
heat is limited. For instance, this scenario would corre- 
spond to post-flare kink oscillations of coronal loops (e.g., 
iNakariakov et al.lll999t [^schwanden et al.l[T999^ . It is un¬ 
likely that such events play a significant role for the heating 
of the solar atmospheric plasma due to the limited amount of 
energy available. However, the results of this paper are also 
applicable to the case of propagating MHD kink waves con- 
tinuously driven by footpoint motions (see, e.g., iTirrv et all 


iDe Ciroof et aklllOOa IPascoe et akllMM l2^ IMTl 
Soleret all 120111) . That scenario would be consistent with 
the ubiquitously observed waves propagating along coronal 
waveg uides (e.g.. lTomczvk et al.l2007trTomczvk & Mclntoshl 
120091) The process of resonant absorption and phase mixing 
operate in the same way for standing and propagating waves 
but, in the case of periodically driven propagating waves, 
there is a continuous energy input and the dissipation of wave 
energy may provide a sustained plasma heating over time. 
This, however, should be checked by using self-consistent 
nonlinear simulations. 

We have compared the actual temporal evolution obtained 
from the modal expansion method with the usual quasi-mode 
approach. In summary, the quasi-mode offers only a partial 
view of the actual behavior of MHD kink waves in nonuni¬ 
form tubes. The quasi-mode provides an effective frequency 
for the global oscillation of the flux tube and a time scale for 
which th e coordinated plasma motion loses its glob al char¬ 
acter (see lRae & Roberta! 198 It ILee & Robertslll986h . but the 
quasi-mode is unable to recover some important features of 
the full temporal evolution. Concerning the apparent damp¬ 
ing of the global kink motion, the quasi-mode is quite accu¬ 
rate when the nonuniform layer is thin, but it becomes less 
and less accurate as the thickness of the layer increases. This 
result may have implications for the accuracy of the seismol¬ 
ogy schemes for kink oscillations that are based on the quasi¬ 
mode, even when the schemes use no approximation and nu¬ 
merically co mpute the exact quas i -mode frequency and damp¬ 
ing rate (see lArregui et al.ll2007t l^ler et al.ll201^ . Another 
limitation is that the quasi-mode does not describe the build¬ 
ing up of small spatial scales and the phase mixing process 
in the nonuniform layer. Therefore, the quasi-mode provides 
no information about how the energy from the global motion 
is redistributed among small spatial scales. The quasi-mode 
becomes a true normal mode in dissipativ e MHD (see, e.g., 
iSteinolfs^l 19851 ITirrv & Goossenslll996h although for real¬ 
istically small diffusivity the dissipative normal mode does 
not have the global character of the ideal quasi-mode and 
is indistinguishable from an ordinyy resistive Alfve n mode 
( jVan Doorsselaere & Poedt^ 120071; ISoler et all 120131) . Hav¬ 
ing these issues in mind is crucial to correctly interpret the 
quasi-mode in physical terms and to reconcile ideal quasi- 
m odes with dissipative eigenmodes (see also the discussion 
in iGoossens et al.ll2014t) . 

In the case of largely nonuniform tubes, the temporal evo¬ 
lution of the transverse displacement of the tube axis reveals 
that, after a first phase dominated by the fast damping of the 
global oscillation, there is a second phase where the oscilla¬ 
tions have a lower but roughly constant amplitude. This sec¬ 
ond phase might be related to the recentl y discovered decay- 
less kink oscillations of co ronal loops (see lNistico et al.ll2013l: 
lAnfinogentoy et ani2013b . Although that relation is specula- 
tiye at present, it might be worth that future works inyestigate 
a possible link between those recent obseryations and the re¬ 
sults discussed here. 
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